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(N ; Abstract 

The search for particle electric dipole moments (EDM) is one of the best places to look for physics 



beyond the Standard Model of electroweak interaction because the size of time reversal violation 
predicted by the Standard Model is incompatible with present ideas concerning the creation of the 
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Baryon-Antibaryon asymmetry. As the sensitivity of these EDM searches increases more subtle 

o' 

systematic effects become important. We develop a general analytical approach to describe a 



systematic effect recently observed in an electric dipole moment experiment using stored particles 



[2]. Our approach is based on the relationship between the systematic frequency shift and the 



velocity autocorrelation function of the resonating particles. Our results, when applied to well- 



known limiting forms of the correlation function, are in good agreement with both the limiting 
cases studied in recent work that employed a numerical/heuristic analysis. Our general approach 
explains some of the surprising results observed in that work and displays the rich behavior of the 
shift for intermediate frequencies, which has not been studied previously. In an appendix we give 
a new derivation of Egelstaf's theorem which we used in our study of the Diffusion theory (low 
frequency) limit of the effect. 
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Introduction 

The search for an electric dipole moment (EDM) of the neutron is perhaps unique in 
modern physics in that experimental work on this subject has been going on more or less 
continuously for over 50 years. In that period the experimental sensitivity has increased 
by more than a factor of 10 6 without an EDM ever being observed. The reason for this 
apparently obsessive behavior by a small group of dedicated physicists is that the observation 
of a non-zero neutron EDM would be evidence of time reversal violation and for physics 
beyond the so-called Standard Model of electroweak interactions. An essential point is that 
the Standard Model predictions of the magnitude of time reversal violation are inconsistent 
with our ideas of the formation of the universe; namely the production of the presently 
observed matter - anti-matter asymmetry requires time reversal violation many orders of 
magnitude greater than that predicted by the Standard Model. 

In this type of experiment (null experiment) the control of systematic errors is of great 
significance. While the switch of experimental technique from beam experiments to ex- 
periments using stored Ultra-cold neutrons (UCN) has eliminated many of the sources of 
systematic error associated with the beam technique, the gain in sensitivity brought by the 
new UCN technique means that the experiments are sensitive to a new range of systematic 



3 



errors. One of the most serious of these is associated with the interaction of gradients of the 
ever-present constant magnetic field with the well known motional magnetic field ( — x E ) . 



As the particles move in the apparatus, these fields, as seen by the particles, will be time 
dependent. This effect was first pointed out by Commins [1] and explained in terms of the 
geometrical phase concept. A more general description is in terms of the Bloch-Siegert shift 
of magnetic resonance frequencies due to the time dependent fields mentioned above [2], [3]. 

The effect was apparently empirically identified in the ILL Hg comagnetometer EDM 
experiment and recently Pendlebury et al [2] have given a very detailed discussion of it, 
including intuitive models and analytical calculations for certain cases, the relation between 
and regions of applicability of the geometric phase and Bloch-Siegert models, numerical 
simulations and experimental verification of the most significant features. However this 
pioneering work has left certain questions unanswered. In particular the understanding of 
effects of collisions on the systematic frequency shifts remains incomplete. 

In this work we attempt to clarify several points concerning the influence of particle 
collisions. We explain the reason that in contrast to gas collisions, collisions with the walls 
were observed to have no effect on the magnitude of the systematic frequency shifts and show 
that this only applies to the limiting cases of high and low frequency. We show that the 
frequency shift is related to the velocity autocorrelation function of the resonating particles. 
Our solution, when applied to well known limiting forms of the correlation function, gives 
results in agreement with those obtained numerically in [2]. McGregor has taken a similar 
approach to the problem of relaxation due to static field gradients [4], whereas the approach 
taken by Cates et al to the problem of static field gradients [5] and gradients combined with 
oscillating perturbing fields [6] is somewhat different than ours. 

A. Brief description of the effect 

Consider a case where, in a storage experiment, there is a radial magnetic field due to a 
magnetic field gradient in the z direction ( B , the quantization axis, and the electric field 
E are along z). Now consider roughly circular orbits, due to specular reflection around the 
bottle at a constant angle, in the x — y plane with radius approximately the bottle radius 
R. The wall collisions occur at a frequency l/r c while the orbital frequency is uo r = 2a /r c 
where a is the incidence angle relative to the surface. We can transform into a rotating 
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frame at u r (note that this is not the Schwinger rotating frame that eliminates -Bo) so that 
the problem is quasi-static [3]. 

The radial field, with the barrel gradient plus v x E field, is 

UOrRE 

B R = B r ±B E = aR± — 

c 

where B r (r) = (r/2) dB z /dz = ar is the radial field due to the axial gradient, and ±B E = 
ru r E/c is the radially directed v x E field and the ± refer to the rotation direction. 
In the rotating frame, 

B 2 = {B -u r /j) 2 + (B R ) 2 

where 7 is the gyromagnetic ratio. Expanding in the limit where Br « B with transfor- 
mation back to the lab frame we find 

l(aR- uj r RE/cf (aR 2 uj r E/c) 
B = B + - = B 



2 B - uj r />y " " B - uj r />y 

keeping only terms linear in Be- Averaging over rotation direction (e.g., the sign of u> r , the 
net effect of the gradient field combined with a v x E yields a systematic (magnetic field) 
shift of 

fa = TSB = (1) 
equivalent to Eq. (18) of [2]. Taking the limit uj r /l <^ B we have 

"--^ *>-- a -W 

where which would seem to set the scale of the effect and is equivalent to Eq. (19) of [2]. 
In this limit, the frequency shift does not depend on 7, implying that it is the result of a 
geometric effect. 

In the other limit, where the rotation frequency is much faster than the Larmor frequency, 
we similarly find that 

5B = ~faR 2 E/c, 5cu = -f 2 aR 2 E/c (3) 

which is independent of the motional frequency u r of opposite sign from the previous limit 
and equivalent to Eq. (21) of [2]. 



I. FREQUENCY SHIFT DUE TO FLUCTUATING FIELDS IN THE X-Y PLANE 
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A. Density matrix approach to the problem 

The issues of the effects of a weak fluctuating potential on the evolution of the density 
matrix have been well-addressed in the literature. However, these treatments generally 
assume that the perturbing potential has a short correlation time, and certain assumptions 
regarding averaging are not applicable to our problem. The effect of a static electric field 
E by itself was treated in [7] where the E 2 effect was related to the correlation time, and 
requirements on the field reversal accuracy were discussed. 

So we therefore start from the beginning, following [8] (p. 276). 

The radial gradient and v x E fields can be treated as weak fluctuating perturbing fields 
B Xj y(t) in the x — y plane, with a constant B applied along z. The perturbing fields B' (t) 
can be written as 

B' x (t) = B x (t) - (B x (t)); B' y (t) = B y {t) + (B y (t)) (4) 

where (....) represents a time average of B x ^ y (t). The constant component of the perturbing 
field are added to B , 

B' = ^j{B Q f + (B x (t))* + (B y (t))* (5) 
leaving the perturbing fields with averages of zero. We define 

^ = 7^; uj x , y (t)= 1 B' x Jt). (6) 

The Hamiltonian is thus 

H = ~^a z - ^a x - ^a y = H + Hl (t). (7) 

Defining 

2b = uj x + iu y ; 2b* = uj x — iujy (8) 
the perturbing Hamiltonian can be rewritten as 

H 1 (t) = b*<r+ + ba- (9) 

where a± are defined in the appendix, and it is understood that b is intrinsically time- 
dependent. Furthermore, the density matrix can be expanded in the spherical Pauli basis, 

P = 1 + Pl,0<*z + Pi,i(*+ + Pi (10) 
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where pn = p\_\- 

The time evolution of the density matrix is 

= -*[#„ + ffi(f),A>]. (11) 

The explicit dependence on the constant H can be eliminated by transforming to the 
rotating frame (also called the interaction representation), with 

H^t) -> e iHo *ifi(t)e- iHo *; p -> e iHot pe- iHot (12) 

where 

( p -iuot/2 n \ 
(13) 
e <«ot/2 / 

We henceforth will work in the rotating frame, with 

H^t) = e~ iu}ot b*a + + e^°V_. (14) 
The time evolution of the density matrix in the rotating frame is 

% = -i[H 1 (t),p] (15) 
which can be integrated by successive approximations to 

p(t)=p(0)-t f[H x {t),p$)W 
Jo 

- fdt' f df[H x (l!),[H x (f),fW)\\- (16) 
Jo Jo 



We are interested in the relaxation rates and frequency shifts due to the perturbing fields, 
which can be found through the time derivative of p, which by introducing a new variable 
r = t- t", is 



dp 
dt 



= -i[ Hl (t),p(0)\- fdriH^^H^t-r)^^)]]. (17) 
Jo 

The first term on the r.h.s has an ensemble average of zero; furthermore, there is no correla- 
tion between p and the fluctuating Hamiltonian (e.g., phases of the neutrons have no explicit 
spatial dependence, and Hi(t) is different for every neutron in the system). In addition, if 
we assume the perturbation is weak, p(0) can be replaced by p(t) which introduces errors 
below second order. 
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We then have 

f t =~ jf dr[#i(0, [Hi(t - r),p(t)}} = Fp(t) (18) 

where T is the "relaxation matrix", the real parts of which describe decays of coherence, 
and the imaginary parts of the off-diagonal elements describe frequency shifts. 

Using the relations in Appendix A together with the expansion of the density matrix Eq. 
(10), the time-derivative of p, correct to second-order and neglecting 2uo terms, is 

= f 2e^b*b'dr (19) 
Jo 

pi,i = -Pi,i [ 2e-^ T bb'*dr (20) 
Jo 

pi,o = -pi,o [ 4Re [e iu)ot b'*b] dr (21) 
Jo 

where 

2b = u x (t) + iu v (t); 2b' = u x (t - t) + iuj y {t - r). (22) 

These equations describe both frequency shifts and relaxations of the density matrix. We 
are at present most interested in frequency shift, which is given by the difference in the 
off-diagonal components of Y. Expanding b and b' we find 

1 /"* 

Su(t) = -- / [cosu t (u x (t)u y (t - t) - u x (t - r)u! y (t)) 

2 Jo 

+ sincj r (uj x {t)uj x {t — r) + ui y {t)ui y {t — r)) ~\dr. (23) 

This is the general solution for the frequency shift given an arbitrary perturbing field. An 
ensemble average must be taken. 

The identical result is obtained with appropriate {w Xty , 5uj uj ) approximations from 
the Bloch equation in the form given in Eqs. (46) and (47) of [2]. This is quite interesting 
given the different assumptions made in the two approaches. 

Now uj x = ax + bv y , ujy = ay — bv x where 



-If 

b = 7- (25) 

c 

with 7 the gyromagnetic ratio and it is clear that only the cross-terms u x u y will result in a 
non-zero linear E (a b) shift, 



8 



1 /•* 

6u = -- dr (cos u„t) {{u x (t) u y (t - t)) - {u x (t -T)u y (t))} (26) 

ah r l 

= ^J0 



dr (cos uj t) R (t) 

^ Jo 
where 

i?(r) = (y(f)u„(* - r) + x(t)v x (t - r) - y(t - r)v y (t) - x(t - r)v x (t)} (27) 
is the net correlation function, where (...) represents an ensemble and time average. 

B. General solution for a radial magnetic field plus vxE 

According to (26) the frequency shift is proportional to the Fourier transform of the 
correlation function R (r) , between (y, v y ) and (x, v x ) evaluated at the Larmor frequency, 
ujq. However this can be written in terms of the velocity autocorrelation function as follows: 

y(t) = vo + / v y (?) dt' 
Jo 

rt—T 

y(t -r)=y + v y (?) d? (28) 
Jo 

Since there are no correlations between y Q and v y the y terms in (27) are 

A = y(t)v y (t-r)= f (vy(?)vy(t-r))d? (29) 
Jo 

B = y(t- r)vy(t) = [ (v y (?) v y (t)) d? (30) 
Jo 



R y (r) = A- B = 



t-T 



= Qf (v v (0 v y (t - r)> d? - ^ (vy (?) v y (t)) dt'j (ab) (31) 
= f dxip (x) - j dxip (x) (32) 

Jr-t Jt 

where ip (x) is the velocity autocorrelation function and we used the fact that it is an even 
function of x. Repeating the same argument for the x axis we have 

^ ( r ) = ( v v (t) v v(t -i~)+v x (t) v x (t - t)) 

= (Vxy (t) ■ ~rfxy(t ~ t)) 
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R(t)= I dxip (x) - I dxip (x) (33) 

Jr-t Jt 

= 2h(r) -h(t-r) -h (t) (34) 
= 2h (r) (35) 

h ( r ) = / dxip (x) (36) 
Jo 

and we consider only cases where ip (x) — > as x — > oo so that we can take the limit 
i — > oo in Eq. (34) and we note that a constant term in R will not have any effect on (26) 
contributing only a term oc 5 (u ) = 0. 

According to (26) we need the cosine Fourier transform of R(t). This will involve Xjio 
times the FT of ip (x) which in turn is proportional to u 2 times the FT of the position 
correlation function as we shall see. Substituting (34) into (26) we have 

5u = ab / dr (cosc<j t) h (r) (37) 
Jo 

Writing the velocity correlation function as 

/oo 
cos utijj (ou) du (38) 
-oo 

we have 

/ sin uit — 1 



o 



h(r)= I i/j(t)dt= I i/j(u>) I ) du (39) 



so that according to (37) the frequency shift is given by (dropping the time independent 
term) 



Suj = ab 



f drcosu o rJ^{u)^dw 



5u = -ab f°° j^^du (40) 

J— OO (U>a—U>' ! ) v ' 

The equation (40) represents the general solution to our problem which is simply the single 
frequency B-S result (Eq. 1, [2] Eq. (18)) summed over the frequency spectrum of the 
velocity autocorrelation function plus oscillating terms (omitted) that don't contribute as 
long as ip (x) — > as x — > oo. 
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1. Example: Particle in circular orbit 

For a particle in an hypothetical circular orbit with orbital frequency u r ^ oj we have 



in agreement with (1 and Eq. (18) of [2]) and valid for the case when (uj — uj r ) > 0J x>y . 
II. NUMERICAL CALCULATIONS OF THE FREQUENCY SHIFT 
A. Numerical estimations of the correlation function 

The problem of the neutron EDM experiment with a 199 Hg comagnetometer subject to a 
time-varying v x E field in combination with a spatially-varying magnetic field is described 
in [2] and in the Introduction. We assume a cylindrical volume with radial field B (r) = a'rf. 
The electric field is constant everywhere and along the z direction. Assuming a constant 
velocity v, the v x E field is then fluctuating in direction but of spatially uniform magnitude. 

A numerical calculation of the correlation function was performed for the two-dimensional 
case (UCN or Hg at a fixed z, moving only in x—y plane). This problem can be parameterized 
in terms of the time between collisions r c = X/v, where the mean free path between collisions 
is A and the average velocity is v. For the numerical calculations, v is assumed constant. 
Time can be parameterized in dimensionless units, r/r c . The correlation function was 
calculated by statistically choosing a propagation distance for a fixed velocity direction, 
and taking time steps of 0.025, after which a new random velocity direction was chosen. 
Various degrees of specularity, parameterized by A9 for the statistical degree of angular 
change for reflection from the bottle surface, were considered. 

Results of a two-dimensional Monte Carlo calculation are shown in Figure 1. Taking 
A = 1 and fixed, we see the effect of wall collisions as the bottle radius approaches A. 



tp( T ) = vl y COS UJ r T 
V'H = vlJ (u - UJ r ) 



(41) 



and substituting in (40) 




(42) 
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We see in Fig. 1 that in all cases R(r) initially increases linearly. The effect of the wall 
collisions when R > A is to limit the distance that the random walk can take, and this 
appears as an exponential decay in R (r) at long times. This effect does not depend on 
the specularity of the wall collisions and is best seen as an effect on the whole ensemble of 
particles which can be described by classical diffusion theory. In this limit, the correlation 
function is well-described by 

R( r ) = (1 - e~ T l Tc )e- T l T (43) 

where, from analysis of the plots, 

T»°4*. (44) 
Xv 

In the other limit, R < A, R(t) oscillates with frequency 

W " MR (45) 

and 

R(r) = e - r/T sin ujt (46) 

where T depends on A0, but is typically of order 27t/lj. 

The frequency shift is determined by Eq. (26) and in the case of large R we find (r c << T), 
using Eq. (44) 

abR 2 abR 2 
= T 2 u 2 + 1 = l + (0.6i?V)/W 
These results are in good agreement with [2], Fig. 10, for which 4/27T 0.634 replaces the 
factor 0.6 above and with Eq. (72) below. 

Additional insight can be gained by considering the effects of varying A keeping R fixed, 
as shown in Fig. 2 for very small A. In this limit, the horizontal axis is multiplied by X/R 
to define time proportional to R/v. The correlation amplitude function is proportional to 
Xv and the decaying exponential time constant is 

r K £. (48, 

The time to reach the peak value is 

t oc X 2 /Rv (49) 

which approaches zero as R — > oo. 
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This limit is further discussed in Sec. 4.1, and the frequency shift in this case is in general 
agreement with Fig. 10 of [2]. 

The curves for large (relative to R) A in figure 1 show damped oscillations whose damping 
depends on the angular spread of the wall collisions. This is a manifestation of the resonance 
behavior discussed in ref. [2] for the case of perfectly specular wall collisions. Here we see 
the damping due to non-specular reflections. 

B. Numerical estimations of the frequency shift for all values of L0 o /^r 

Using Eq. (37), and the results of the previous section, the cosine transform of the 
numerically-determined correlation function can be calculated numerically. In order to re- 
duce oscillations due to the finite time window, a Hamming window function was applied 
to the correlation function, and a slight correction due to the frequency dependent gain as 
imposed by the window function was applied. The results, as a function of mean free path 
A at fixed radius R, for specular and purely diffuse wall reflection, are shown in Fig. 3. 

There are a few points worth noting. First, the curves for large A in the specular case are 
very similar to the Bloch-Siegert result. Second, at small and large frequencies, the results 
agree with the numerical semi-analytically determined results presented above, and in [2] 
and the theoretical analysis below. Third, the behavior at intermediate frequencies is seen 
to be very interesting: The shift goes to zero for 0J o /u r ~ 1 as it must because the effect 
changes sign between large and small frequencies. 

Furthermore, it can be seen immediately that the effects of wall collision specularity 
is important when u Q ~ u r , in contradiction to the statement in [2] that the degree of 
specularity does not affect the frequency shift. We discuss this point later in more detail 
(Sec. IV). 

III. ANALYTICAL RESULTS FOR THE LIMITING CASES OF LARGE AND 
SMALL FREQUENCIES (u /u r > l,u a /u r < 1) 

Equation (40 ) represents the formal solution of the problem in all cases of interest here. 
Thus the frequency shift is determined entirely by the velocity auto-correlation function of 
the particles undergoing magnetic resonance. This function has been the subject of intense 
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experimental and theoretical study ([13-15]). In our case, involving macroscopic distances 
and times, it suffices to treat the motion classically. For relatively short times if the particles 
undergo collisions which are distributed according to a Poisson distribution with average time 
between collisions given by r c , the velocity correlation function is well known to be given by 

^ (t) = (v*) e -*/rc (50) 

This form is known to be valid for relatively short times. According to Eq. (37) the 
frequency shift depends on the Fourier transform of the integral of the velocity correlation 
function evaluated at u . So the short time behavior of ip (t) determines the high frequency 
behavior of ip (cu), and the result using this form is expected be valid in the case of large u a , 
i.e. ijj ^> uj r . 

For longer times the velocity correlation function is well described by classical diffusion 
theory. Thus the long time behavior will determine the low frequency region of the velocity 
spectrum and the result will apply to the case u Q <C uj r . In this region the result will depend 
on the size of the containing vessel as the dynamics of the diffusion process are influenced 
by the boundary conditions. 



A. Short correlation times (io r <C lu ) 



Using (50) we have 

1 r°° 1 1 

^(u)^-(v 2 ) COS ute-^ T dt= -(v 2 ) -. zr-\ (51) 

7T J 7T T {CO 2 + 4j) 



so that according to (40) 



5uj = -ab I , f( U \. du (52) 

k 2 - ^ 2 ) 



oo 

2\ r°° 



1 (v 2 ) 

^—L ^+Mw (53) 

_ a6 M^_L^ (54) 



< 1 + J 



2 T 2 



This is in substantial agreement with the expression given in the caption of [2] Fig. 12 when 
it is taken with [2] Eq. (19) or (3) applicable to the case when u r <^ uj . It is quite likely 
that the small discrepancy (~ 10%) in the 50% suppression point is due to the process of 
averaging over the velocity distribution in [2] Fig. 12. 
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B. Diffusion theory calculation of the long time behavior of the velocity correlation 
function. Frequency shifts for (uj r » u> ) 

Whereas the previous case applies to UCN this case would apply to atoms used as a 
comagnetometer and is more relevant experimentally as it results in larger shifts [2] and 
in some cases [9] the collision rate can be simply adjusted by changing the experimental 



In the following we review the solution of the diffusion equation in cylindrical geometry, 
obtain the velocity autocorrelation function from the solution and calculate the frequency 
shift. In the limit of small collision rate the result agrees with the known results for {uj r ^> uj ) 
and the effect of the collisions agrees with that found from numerical simulations ([2] Fig. 



1. Green's function for the diffusion equation in cylindrical geometry 

In this section we attempt to understand the effects of the vessel boundary on the velocity 
autocorrelation function, observed in the numerical simulations (section II A), by applying 
classical diffusion theory to the problem. Diffusion theory is expected to be valid for long 
times so that we expect the results to be valid for small cu , i.e. cu <C uj t . 



We consider a two dimensional problem, that is we neglect any z dependence ik z = 0) . 
For the cases considered in [2] where the height of the bottle is much smaller than the radius 
higher z modes will decay relatively quickly. 

The boundary condition is j (R) = — Dt£ = so the eigenfunctions satisfying the bound- 
ary conditions are 



conditions. 



10) 





V 2 m + k 2 u = 



(55) 




k m , n R = x' mn ( nth zero °f dJ m (z)/dz) 
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where the normalization constant (which depends on the boundary conditions) is ([10], p 
322) 

N 1 / 2fc m,n 

\[ 7 hsJ m {k m , n R) y {k m . n R) 2 - m 2 
The Greens' function satisfying the boundary conditions is [16] 

G(r, r', ( N ™,n) 2 JUk^JUk^y^ 9 - 9 ^- ^ (57) 

m,n 

This is the probability of finding a particle at ~f at time t, given that the particle was at 
~r' at time t — 0. The spectrum of the velocity correlation function is related to S(~q,uj) 
which in turn is the average over the system of the Fourier transform of this probability with 
respect to p = (r — r'). We use the cosine transform because we want the cosine transform 
of the velocity correlation function (38) ip (cu). 



(58) 



S(q, uj)=^-lj d 2 pS^ dt cos utG(r, r', t) 

Iff d 2 r' -> -> f°° 
= - J J — J dtcoswtG(r,r',t) (59) 

= ^R2 E ( N m,nf J d 2 re^ J m (k m , n ry me x (60) 

m,n 

/Dk 2 
d>r'e-^'j m (k m>n r>)e-^' (61) 



Dk 2 

Now we can evaluate the integrals using 



J m ( x ) = i-^/ e i(*ooe«+m«) d0 (62) 
27T 





and Bessel function identities 



/c? re 9 Jm(k mn r^e , . Jm(k m ^ n R)x 

(Q 2 ~ k 2 ) 2 



(<Z 2 -^,n) 2 
[J m _i(gi2) - J m+ i(gi?)] 



thus 

^) = -f^E" (63) 



m,n 



k'kn ( 2tt gi? \ 2 -DA: 

[J m _i(gi?) - J m+ i(gi?)J 



rn.n 



((k m , n R) 2 - m 2 ) \(q 2 -kU 2 L ™ ' m ™ ^ (W(^U 2 ) 

(64) 
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2. Velocity autocorrelation function 



The velocity autocorrelation function 

iP( T ) = ( "?/ (0) • "7/ (r)) (65) 
has a Fourier transform given by ([11]) 

^H = hm2^ S(q,u>) (66) 

so that the only terms in (64) which contribute are those containing J (qR), since 
lim^o Jn (x) ~ (x) n , J (0) = 1. Thus we only need to keep terms with m — ±1 in (64) and 
we find 

\\mS(q,Lo) — — ^5 -—, r- (67) 

n ((h, n R) 2 - 1) (u,2 + ( Dkln) ^ 



Then 



^H = -y /9 4 , 7 — — — ^ (68) 



5. Frequency shift in the diffusion approximation (cylindrical geometry) 
According to (40) 



1 

5uj = —ab I i])(uS)— -du (69) 



u 2 - u 2 



= abR 2 ^2 7 2 ^ — 77 7 " — 

n [ X l,n - 1 ) 2 / ( avgi \ 



2 

+ 1 



(70) 



The result (70) is dominated by the first mode x 1:1 = 1.84. Figure 3 shows the first term 
in comparison to the sum of the first 4 terms. For convenience we list the zeroes of J[(x): 
xi t 2 = 5.33,2:1,3 = 8.54,2:1,4 = 11.7. Since we are dealing with a 2 dimensional problem we 
put 

D = v 2 t/2 (71) 
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(instead of tv 2 /3 for 3 dimensions) in order to facilitate the comparison with the numerical 
simulations [2] and obtain for the condition that the frequency shift is reduced to 50% of its 
value in the absence of collisions 



uj R 2 2u D R 2 u R 2 

1] = - a = — — = -59 = 1 



V 2 T 



(72) 



Dx\ n v 2 Tx\ n 

the numerical factor of which is to be compared with | = .634 obtained in [[2]], fig. 10 
by fitting simulated results, and our numerical result of 0.6 presented in Sec. II A. The 
magnitude abR 2 /2 of (70) in the absence of collisions is just that expected from the Bloch- 
Siegert treatment in the case u r 3> uo (Eq. (3), [2] Eq. (21)), averaged over the different 
trajectories as discussed in [2] after equation (22). 



4- Frequency shift in the diffusion approximation (rectangular geometry) 
For the rectangular case the normalized eigenfunctions are 




L T l 




2 ran / 2 mir 
cos — — X\ I — cos y 



Li. 



(73) 



J y 

which satisfy the reflection boundary conditions at x — 0, L x and y — 0, L y . For n or m = 
the corresponding eigenfunctions are 

1 



u = 



(74) 



x,y 



so that the Green's function is 



G(x,x',y,y',t) = ^ 



m,n=0 



12 , _£)^2 t 

~ ~~\~~ ~ COS KjyiX COS K m X 6 



x 



1 2 

- — h — cos k n x cos k n x e 

Ly Ly 



(75) 



with k mjn = (m, n) ir/ L x y . To calculate lim^o S (q, u) we need integrals of the form 

lim / e iqxX cos k m xdx 
Jo 

= {L X m = | 



q x L x m = 2,4,6.. 
" k 'i \ =1 m = 1,3,5.. 



(76) 
(77) 



Since each of these will appear squared because of the contribution from the x, x' integrals 
we can only take the odd values of m. The even numbers will yield a higher power of q 
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which will vanish in the limit. Given this, if we take m — 1, 3, 5 we must take n = and 
vice-versa. We calculate, using (76) 



lim S(q,w)= fdx f^- fdy f ^ e ^-(W)I f° dt coscotG (x, t' , t) 

1^ l J J L X J J Ly 71 J 



Dkl 



Em=l,3, 

^"=1,3,5.. ki(nwf ui 2 +(Dklf 



Dkl 



where we used (ql) = (q^) = q 2 /2. Then 



E 



Dkl 



+ E L.2 („, 



Dki 



and (using 40) 



f°° 1 
5u = —ah I -duo 

( 



= 8ab 



E 

m=l,3,5. 



1 



[mix 



T 2 

+ E " 



\ 



yD(m7r) / 



+ 1 



/ 



(78) 
(79) 



(80) 



(81) 
(82) 



We thus see that in a rectangular box L x ^ L y it is the longer side which dominates the 
behavior. 



C. Diffusion theory calculation of the long time behavior of the velocity correlation 
function. Frequency shifts for (oj r » 0J o ) (Alternate caclulation) 

In this section we derive the diffusion theory result for cylindrical geometry (70) using an 
alternate method based on that of Mcgregor [4] which avoids the use of the theorem (66). 
We start with the Green's function for cylindrical geometry given above (57) 

G(r, r', f)=£ {N m , n f J m {k m , n r) J m {k m>n r')e im ^ e~ Dk ^ (83) 

m,n 

This is the probability of finding a particle at ~f at time t, given that the particle was at 
~r' at time t — 0. 
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1. Position-velocity correlation function R{t) 



From equ. (27) we have 

R(r) = (y(t)v y (t - r) + x(t)v x (t - r) - y(t - T )v y (t) - x(t - r)v x (t)) (84) 
= ( "r (t) • 1?(f - r) - ~?(t - r) • 



^ (?(t) ■ ?(tr)) - j t (?(t/) ■ ?(t)) 



t'=t-T 

(85) 



Following McGregor [4] we write 



(?(t) ■ r (*/)> = J d 2 r j ^G(r, r>, t - t') (? ■ 7) (86) 
= ^Y1 ( Nm ^ J rdrd(f>J m (k m , n r) J r'dr'd<j)'J m (k m , n r') x (87) 

m,n 

x rr'e^-^ [cos cos 0' + sin sin 0'] e - Dfe ™.«(*-*') (88) 

Since (N mtn ) 2 = (iV_„ ti „) 2 and J- m (z) = (—l) m J m (z), (56) we can (excluding m — 0, which 
will be seen not contribute to the result) combine the terms for m and — m as follows: 



(7>(f) • 7>(t/)) = Yl ( iVm '«) 2 / r2drJ m(km, n r) [ r' 2 dr'J m (k m>n r')x (89) 
* J d<p J rf0'S(0,0 , )e- Dfc -." ( '-'' ) (90) 



where 



S(0, 0') = [cos cos 0' + sin sin 0'] 2 cosm (0 — 0') (91) 

= 2 [cos cos 0' + sin sin 0'] [cos m0 cos m0' + sin m0 sin m0'] (92) 

= - (cos cos m(f) cos 0' cos m(f)' + sin sin m0 sin 0' sin m0'+ (93) 

+ sin cos m(f) sin 0' cos m0' + cos sin m0 cos <f>' sin m0') (94) 

= - [(cos (m — 1) + cos (m + 1) 0) (cos (m — 1) 0' + cos (m + 1) 0') + (95) 

+ (cos (m — 1) — cos (m + 1) 0) (cos (m — 1) 0' — cos (m + 1) 0') + (96) 

+ (sin (m - 1) - sin (m + 1) 0) (sin (m - 1) 0' - sin (m + 1) 0') + (97) 

+ (sin (m — 1) + sin (m + 1) 0) (sin (m — 1) 0' + sin (m + 1) 0') (98) 
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using 

»2tt r2ir 



/ d(j)cos(m + 1)0 = (m>0), / d0sin(m ± 1)0 = (99) 

Jo Jo 

/ d0cos(m - 1)0 = 27r5 ml (100) 
Jo 



we have 

\2 



j d<p j d0'~(0,0') = ^-25 ml (101) 



so that 



(?(t) ■ ?(tt)) = ^ (2vr) 2 £ 5 ml (N m>n ) 2 J r 2 drJ m (k m , n r) j r' 2 dr'J m (k m , n r')e- Dk ^-^ 



Using 



we have 



Thus 



and 



m>0,n 

(102) 



= ^ E (^i.™) 2 £ r-'dr-Uk^r) j\' 2 dr' J^h^e^U^') 



n 



(103) 



I" z n dzJ n ^(z) = \z n J n (z)]l (104) 
Jo 



= -J n (z) - J n+ i(z) (105) 

z 

f R r'drMh^r) = J^J^J (106) 

<"?(*) • r>(*)> = ^E (^^) 2 7 ^4^(^,n)e- Dfcf - (t - t ' ) (107) 

= 4jR 2 y 1 e -^f,J*-*') (log) 

rKJ a (KJ a -i) 



R(r) = 2j2vnDkl n e- Dk l^ (110) 
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using (85). Then according to (26) 




dr (cos uj t) R (t) 



(111) 



ab / dr (cos u t) rj n Dkl n e Dk ^ T 



(112) 




+ 1 



(113) 



in complete agreement with equation (70). 
D. Application: 3 He Comagnetometer 

In [9] the use of 3 He as a comagnetometer for a UCN neutron EDM experiment is dis- 
cussed. This system is rather unique in that an effective background gas (phonons) can be 
introduced which affects the 3 He significantly while having no substantial interaction with 
the UCN for temperatures below 0.5 K. Because the 3 He and neutron magnetic moments 
are equal to within 10%, it is possible to control this systematic by varying the size of the 
effect for 3 He by changing the diffusion rate of the 3 He. 

The UCN upscattering lifetime varies as 100T -7 s for T < 0.7 K, while the coefficient of 
diffusion for 3 He in a superfluid helium bath varies as D ~ 1.6T~ 7 cm 2 /s [17]. 

In connection with (72) this yields rj — 1 when the superfluid helium temperature is 
T m 0.25 K, (R=25cm), which determines the temperature scale where the effect can be 
varied, and is within the design range of operating temperature for the planned experiment, 
compatible with a UCN upscattering lifetime in excess of 1000 s. 

IV. DISCUSSION 

One of the surprising, but unexplained results of [2] was that according to their numerical 
simulations, wall collisions had no influence on the magnitude of the frequency shifts while 
gas collisions could eliminate the frequency shifts completely if their rate is high enough. 
This was apparently only studied in the limits of large and small u /uj r . We now know that 
this does not apply to intermediate frequencies, e.g, when u ~ u r . In Fig. 3 we see that 
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wall collisions have a serious influence at intermediate frequencies when A > R. Also from 
Fig. 3 we see that the curves for diffuse wall reflections in the absence of gas collisions is 
very similar to the specular curves for A < R/2. This implies that there is no essential 
difference between wall and gas collisions. We now show that the reason the wall collisions 
have no effect at the limiting frequencies, contrary to the case at intermediate frequencies, 
is that the wall collisions are never fast enough to influence the systematic (proportional to 
E) frequency shifts in the limits of large and small uj . 

For a particles in a cylindrical vessel following a trajectory along a chord subtending an 
angle 2a, the time between collisions is 

2R . . 

t c = — sin a (114) 

and the effective field rotation frequency is given by 

av 

u r = 2a t c = — ; (115) 

it sin a 

Considering first the case when u r ^> u (70, [2] Fig. 10,) the systematic frequency shift was 
found to be suppressed by the factor r\ 

1 

rj 



l + (3 2 
2i?V, 

P n 

TTV Z T C 

For significant suppression we need (5 > 1 

2R 2 uj Ruo 



7TV 2 T r irv sin a 



> 1 (H6) 



Rw _ 25 x 2 • 7 1 

^ ~ lb 1 30 

for representative conditions in [2], fig. 10. (R = 25, B Q = lfiT,v = 10 4 cm/ sec). 
The probability of a given value of a is given in Eq. (Bl) of [2] as 

4 

P(a)da = — sin 2 ada (117) 

7T 

P(a < e) ~ e 3 

so that the wall collisions would only be expected to be effective for a vanishingly small 
fraction of the trajectories. 
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Turning now to the case uj r <C co (54, Fig. 12 of [2]) we have as the condition that the 
suppression be effective: 



1 . 2R 



sin a < 



> — sin a 

L0 o V 

200 



2Ru 2 • 25 • 200 50 

for conditions typical of [2] Fig. 12 (v = 200cm/ sec,B D = 1/jlT). 

Thus the wall collisions rate is never high enough to significantly effect the magnitude 
of the frequency shift at the limits. The wall collisions do, however, broaden and shift the 
resonances discussed in [2] 



V. CONCLUSION 

We have developed a general technique of analyzing the systematic effects due to a com- 
bination of an electric field and magnetic gradients as encountered in EDM experiments 
that employ gasses of stored particles. Use of the correlation technique, either by numerical 
calculations for complicated geometries, or by the velocity correlation function for simpler 
geometries, provides a simplified approach to the problem compared to numerical integra- 
tion of the Bloch equations. Our analysis has added insight to this new systematic effect 
and provides a means of rapidly assessing the effects of various geometries and angular 
distributions for wall and gas collisions. 
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VII. APPENDIX 1: MATRIX ALGEBRA OF SPHERICAL PAULI MATRICES 



The following relationships among the Pauli matrices have been employed in the calcu- 
lation in section I A. 

2<7-|- = a x ± iijy (118) 
o±o" 2 = Tcr±; a z a± = a ± (119) 

v±°t = \ ± \ a * ( 12 °) 
a z a z = 1; o±o± = (121) 



VIII. APPENDIX 2: EGELSTAFF'S VELOCITY CORRELATION FUNCTION 
THEOREM; A NEW LOOK AT AN OLD THEOREM 

A. Introduction 

The relation between the velocity autocorrelation function (vacf) and S s (q,ou) which we 
used in section [III B] was first introduced by Egelstaff [18] and has proven to be a useful 
tool in the study of liquids. The vacf can be simulated for various models and obtained 
from neutron scattering data using Egelstaff 's theorem. The theorem has been discussed by 
several authors [13], [14], [15] and has been given in slightly different forms depending on 
the normalization chosen for the functions involved. 

Following Squires' [15] derivation would yield 

^(oo)=3Lo 2 \im Smc{q ' Uj) (122) 

q^O q 2 

if we were to define 

^ M = i/°° ( ^ (0) '^ (r))e "^ r 

which has a different normalization then used by Squires. Egelstaff gives the theorem as 



where he defines 



q 2 



^ f°° (v x (0)-Mr))e- tuJT dr 
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which accounts for the factor of 3 difference. Both authors give the derivation only for the 
Gaussian approximation where we take 

_ 1 - 2 
G(r ,t) = —-e 2 -M 

(27rw(r)) 3/2 

J(g ; r) = J rf 3 rG(T > ,t)e^- r = e' q2 ^ (123) 

Boon and Yip [19] derive the theorem for the general case, i.e. without the Gaussian 
approximation. 

The theorem has been used to extract vacf's from neutron scattering data by many 
authors. An early example is given by [18]. See also the work of Carneiro [20]. 

The fact that sections III B and III C give the same result and section III B uses Egelstaff 's 
theorem while section III C does not, suggests that we have discovered a new way of proving 
the theorem. 



B. A new derivation of Egelstaff 's theorem 

In this section we will give a general derivation (not relying on the Gaussian approxima- 
tion) of Egelstaff 's theorem. 

We begin by following the formulation of Squires [15] and calculate the velocity autocor- 
relation function 

^(r) = (V(0) • l?(r)> 

as follows: 

Let ~r(t) be the position of a particle at time t, when the particle was at the position 
7^(0) at time t — 0. Then 

r(t) 4 ~r (t) - r (0) = f ~v(t')dt' 

Jo 

r 2 (t) = [ l?(t')dt' ■ [ H(t")dt" 
Jo Jo 

= 2 f dt" f {l?{t')-l?{t"))dt t 
Jo Jo 

Now since (~v(t') ■ ~v(t")) = f(t" — t') (for stationary systems) we can write 

r 2 (t) = 2 f (v(0) ■ ~v(t')) (t - t') dt' 
Jo 
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and 



d 
~dt 



r 2 (t) = 2 (1?(0) • ~v(t)} t + 2 f (~u(0) ■ ~v(t')} dt' - 2 (1?(0) • "??(*)) * 

Jo 

= 2 [ (v(0)-l?(t'))dt' 
Jo 

^r\t) = 2 0?(0) • T?(f)> (124) 

Now, based on the usual definition of the pair distribution function, G(~r,t), we have (fol- 
lowing 



(t) = y r 2 G(r,t)rf 3 r 

| d 3 ql{^,t)^ 



.2 ^ / 



(27T) 3 

= | d 3 ql(t,t) J -^r^- 



-J d 3 qI(l?,t)V q -(v q J ^ 



d 3 r 1=> 1. 
(2^f7 

=-/ rfV(t,ov 9 - (^(t)) 

where we introduced the spatial Fourier transform of the pair distribution function, I(~q, t). 
(See the first part of equation 123) 
Integrating by parts twice we obtain 

r 2 (t) = J d 3 qV q I(q,t)-(v q 5 3 (lf)) 

= _y d 3 qS 3 (t) V 9 - (y q I(q,tj) 
= ~ S ' (^V(?> 0) = - Jim (V 9 2 /(g, 0) 
Then, using (124) we have 

1- (2d 2 



^(r) = (7(0) • v(t)) = --Urn f V 2 — I(q,r) ) (125) 
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Writing 

/CO 
e iuJT tP (u>) duo 
-co 
pco 

J(g ?r ) = / e iuJT S(q,tu)dtu 



we find (in general) 



^H = yfm(vl%o;)) (126) 

This appears to be different than the usual form of the theorem (122) but is completely 
equivalent as can be seen by using 



S 8 (q,u)=Jdte- iut (/?-(3(*)-3(0)) 
and expanding for small q. 

f .J _ - x • (ify) " %)) 
S s = / die-*" il + fq- yR(t) - R(0)) - ± ^ ^ 



The first term gives a 5 (a;) which does not contribute to the result, the second term averages 
to zero and the third term gives 

S s (q,u) = aq 2 

(for isotropic media). Then 

Km (v 2 q S(q,u))= 6a 

and equation (126) is equivalent to (122) in general. 
For the Gaussian approximation 

I(q,r) = e q 2 

we have 

d 2 Tl . q 2 d 2 w 
hm — — i (o, r) = — ■ 

q^odr 2 w ' J 2 dr 2 

and using (125) 

1,. ( 2d 2 rt C\ ?>d 2 w ,.3 d 2 Tl , 
^( T) = - 5 bm(v ? -/(,, T )j =-— = - --/(,, r) (127) 

Finally 

3 

V> (w) = lim — cu 2 ^(g, w) (128) 

for the Gaussian approximation, (note that V^g 2 = 6), confirming that our formulation again 
gives the correct result. 
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C. Discussion 



For the case of diffusion, at very long times the Gaussian approximation holds (w(t) = 
2Dt), but d 2 w/dr 2 = 0, so that the calculation in (127) appears to break down. Nonetheless 
the more general derivations show that one can apply Egelstaff 's theorem to this case where 

1 Dq 2 
7iu 2 + (Dq 2 Y 

as was done in section III B. 



S(q,<*) = z o , L ox2 - ( 129 ) 
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X. FIGURES 



2F 




FIG. 1: The position-velocity correlation function, R(t) = 2h(r) as a function of cell radius R 
parameterized in terms of the mean free path A for different degrees of specularity as parameterized 
by the angular spread of the final angle compared to the incident angle. Red: specular; Green: 
45°; Blue: 90°; Black: diffuse. For R > 2.5A, there was practically no effect due to the degree of 
specularity, as expected. 
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z/IR/v] 



FIG. 2: The position-velocity correlation function when A is very small. This represents the limit 
for slow diffusion. 




FIG. 3: Results of numerically applying Eq. (37) to numerical calculations of the correlation 
function, for varying A with R fixed. 
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